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Introduction 


MicroRNAs  (miRNAs)  are  recently  discovered  but  minimally  explored  non-protein- 
coding  double-stranded  RNA  molecules  that  regulate  gene  expression  involved  in  the 
control  of  development,  proliferation,  differentiation,  apoptosis,  stress  response  and 
cancer  (2).  Recent  studies  of  miRNAs  in  numerous  malignancies  have  demonstrated  that 
they  are  frequently  located  at  genomic  regions  involved  in  cancers,  play  a  critical  role  in 
post-transcriptional  regulation  of  oncogenes  and  tumor  suppressor  genes  (4),  and  are 
associated  with  poor  prognosis  in  leukemia,  colorectal  neoplasia,  lung  cancer,  and  breast 
cancer  (1-3,  5).  The  functional  roles  of  miRNAs  in  breast  cancer  etiology  and  tumor 
progression  are  not  understood. 

Relevant  mouse  models  of  human  breast  cancer  have  been  successfully  used  as  a  filter  to 
help  identify  evolutionarily  conserved  gene  expression  patterns  in  breast  cancer  (5).  It  is 
quite  likely  that  genes  and  genetic  networks  that  are  found  to  be  operative  in  both  the 
mouse  and  human  cancers,  are  critical  to  the  development  and  maintenance  of  the  tumor 
phenotype.  Use  of  mouse  models  in  conjunction  with  human  tumor  data  provide  a  means 
to  overcome  inherent  limitations  of  expression  variance  in  the  human  population  due  to 
great  genetic  diversity  in  the  population.  This  approach  can  leads  to  the  identification  of 
high  value,  functionally  conserved  genes.  The  identification  of  such  genes  also  leads  to 
the  identification  of  candidate  genes  for  targeted  therapies. 

We  have  demonstrated  the  relevance  of  several  existing  mouse  mammary  cancer  models 
that  can  be  directly  clustered  with  three  sub-types  of  human  breast  cancer  using  the  Sorlie 
classification  gene  signatures  (Green,  unpublished  data,  and  7).  This  proposal  uses  this 
same  principle  of  integrating  mouse  model  data  to  identify  the  evolutionarily  conserved, 
functionally  most  important  miRNAs  related  to  human  breast  cancer.  By  integrating 
studies  of  relevant  mouse  models  with  human  breast  cancer,  this  approach  should  lead  to 
the  identification  of  the  likely  miRNA  candidates  that  are  functionally  important  for 
breast  cancer  development  and  that  may  serve  as  new  targets  for  breast  cancer  detection, 
diagnosis  and  therapy.  Functional  validation  of  the  candidate  miRNAs  can  be 
subsequently  performed  in  in  vitro  and  in  vivo  assays. 
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Body 

Hypothesis/Rationale/Purpose. 

An  integrated  analysis  of  miRNA  expression  in  GEM  mammary  cancer  models  with  that 
of  human  breast  cancer  will  provide  significant  insights  into  the  roles  of  miRNAs  in 
oncogenesis.  A  cross-comparison  of  human  and  mouse  miRNA  data  will  identify 
evolutionarily  conserved  candidate  oncogenic  miRNAs  that  likely  are  functionally  highly 
significant  for  breast  cancer  pathogenesis.  These  candidates  can  be  functionally  tested  for 
their  role  in  tumorigenesis,  may  serve  as  potential  biomarkers  and  may  identify  new 
targets/pathways  for  therapies. 

In  order  to  more  fully  compare  genomic  data  between  mouse  and  human,  we  have 
additionally  initiated  array  studies  to  assay  a  larger  number  of  genes  (than  present  in  our 
older  datasets)  using  an  Affymetrix  platform  and  array  CGH  to  identify  alterations  in 
copy  number  changes  in  chromosomal  DNA  using  an  Agilent  platform. 


Experimental  Approaches. 

Collection  of  mammary  tumors  from  mouse  models. 

Multiple  tumor  samples  were  collected  from  several  animal  models  including  C3(l)/Tag, 
MMTV-her2/neu,  MMTV-myc,  MMTV-PyMT;  p53  ""ER+  transplant  model,  p53 "  "ER- 
transplant  model,  and  BRCA1  ;  p53H  '  model.  MMTV-ras  tumors  were  not  initially 
available  and  animals  have  been  bred  to  generate  these  tumors  for  inclusion  into  the 
study. 

Isolation  of  miRNA.  mRNA  and  DNA  from  the  tumor  samples. 

Total  RNA  was  isolated  from  the  tumor  samples  using  a  method  that  did  not  exclude 
small,  miRNA  species.  This  allowed  for  the  use  of  the  total  RNA  for  both  miRNA  and 
expression  profile  arrays.  DNA  was  also  extracted  from  the  same  tumor  samples  for 
array  CGH  analyses.  RNA  and  DNA  were  quality  controlled  using  an  Agilent  Lab-on-a- 
Chip  Bioanalyzer.  Samples  with  poor  quality  were  excluded  from  analyses  and  have 
been  or  will  be  replaced  to  fully  populate  the  datasets. 

miRNA  array  analyses 

Total  RNA  from  tumor  samples  were  labeled  with  cy-3.  Ambion  mouse  reference  RNA 
was  used  as  the  standard  reference  RNA  since  it  is  produced  in  a  manner  that  captures  the 
miRNA  species  and  was  labeled  with  cy-5.  Initially,  we  evaluated  several  miRNA 
platforms  to  determine  which  platform  would  produce  the  best  results.  We  found  that  the 
Ambion  miRNA  platform  was  not  as  sensitive  as  other  platfonns  and  therefore  decided  to 
conduct  all  analyses  using  an  Agilent-produced  array  that  was  designed  by  the  Frederick 
miRNA  array  core.  This  platform  was  designed  as  a  2-channel  platform  that  gives  highly 
reproducible  results  with  good  sensitivity  and  specificity.  Results  of  miRNA  analyses  are 
provided  in  Figure  1  and  Figure  2. 
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miRNA.AnovaList103 


Figure  1. 

Differentially  expressed  genes  were  identified  with  one-way  ANOVA  (univariate  F-tests) 
across  all  models  (7-level  factor).  The  genes  differentially  expressed  between  the  models 
were  found  with  the  alpha  level  cutoff  at  0.01.  Hierarchical  clustering:  the  103  genes 
were  median  centered  and  the  1  minus  Pearson  correlation  distance  with  complete 
linkage  were  applied  to  generate  the  dendrogram. 
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Figure  2. 

Hierarchical  clustering:  the  103  genes  as  described  in  Figure  3  were  median  centered  and 
the  1  minus  Pearson  correlation  distance  with  complete  linkage  were  applied  to  generate 
the  heatmap.  The  clustering  shows  that  her2  and  c-myc  have  the  most  defined  miRNA 
expression  profiles  (within  cluster  correlation  around  0.6). 


Expression  profiling  analyses 

All  tumor  samples  were  analyzed  using  the  Affymetrix  U130  array  which  provided 
significantly  more  gene  representation  than  data  previously  generated  in  our  lab,  thus 
providing  additional  data  to  correlate  with  the  miRNA  data.  Since  the  targets  of  miRNAs 
are  poorly  understood  and  have  generally  only  been  identified  by  computational  methods, 
but  not  validated  biologically,  this  additional  level  of  information  will  be  valuable  in 
selecting  miRNA  targets  for  further  validation.  Data  from  these  studies  indicate  that  the 
tumor  models  cluster  into  separate  categories  as  expected  (Figures  3,  4  and  5).  The 
clusters  based  upon  expression  can  be  compared  with  the  clusters  identified  by  miRNA 
array  analyses. 
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Figure  3. 

Differentially  expressed  genes  were  identified  with  one-way  ANOVA  (univariate  F-tests) 
across  all  models  (7-level  factor).  The  genes  differentially  expressed  between  the  models 
were  found  with  the  alpha  level  cutoff  at  0.001.  Hierarchical  clustering:  in  each  case  the 
genes  were  median  centered  and  the  1  minus  Pearson  correlation  distance  with  average 
linkage  were  applied  to  generate  the  dendrograms.  Gene  selection:  4949  probes  to 
generate  sample  dendrogram. 
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Figure  4. 

500  top  genes  (p  value  rank)  to  generate  two-way  clustering  heatmap.  The  hierarchical 
clustering  analysis  shows  clear  model-specific  signatures,  which  are  robustly  preserved 
across  the  gene  selections  described  above. 
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Figure  5. 

500  top  genes  (p  value  rank)  to  generate  two-way  clustering  heatmap. 


Array  CGH  analyses 

DNA  from  all  tumor  samples  has  been  used  to  perform  2-channel  array  CGH  studies 
using  the  high  density  Agilent  mouse  CGH  array.  It  is  clear  that  alterations  in  genome 
copy  number  vary  between  the  models  and  this  may  be  related  to  changes  in  miRNA 
expression  and  gene  expression.  See  appendix  for  representative  information. 

Functional  analyses  and  bioinformatics  analyses 

The  complexity  of  genomic  alterations  in  the  models  studied  have  led  to  significant 
informatics  challenges  that  we  are  currently  addressing.  No  tools  are  readily  available 
that  allow  us  to  cross-compare  expression  profiles  with  miRNA  profiles  and  array  CGH 
data  together.  We  are  working  on  methods  to  improve  these  types  of  analyses.  This  has 
delayed  our  completion  of  identifying  all  of  the  miRNA  candidates  we  will  functionally 
test  for  their  role  in  the  development  and  progression  of  breast  cancer.  This  work 
remains  ongoing  and  a  high  priority. 
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Key  Research  Accomplishments: 

1)  collection  of  tumors  from  multiple  mouse  models  of  human  breast  cancer  initiated 
by  different  oncogenic  events; 

2)  determination  of  miRNA  expression  from  these  models  and  cross-comparison  of 
expression  between  the  models; 

3)  determination  of  gene  expression  profiles  from  the  same  tumors  and  cross¬ 
comparison  of  the  expression  profiles  between  the  models; 

4)  determination  of  genome-wide  copy  number  alterations  using  array  CGH  and 
cross-comparisons  between  the  models; 

Reportable  Outcomes: 

1)  miRNA  expression  patterns  are  quite  distinct  between  different  models  dependent 
upon  the  initiating  oncogenic  event; 

2)  MMTV-her2  and  MMTV-myc  models  appear  to  have  the  least  number  of 
alterations  in  miRNA  expression  patterns; 

3)  The  alterations  in  miRNA  expression  in  the  MMTV-her2  and  MMTV-myc 
models  are  generally  distinct  from  those  observed  in  the  C3(l)/Tag,  MMTV- 
PyMT,  p53-/-  and  BRCAl-/-;p53+/-  models; 

4)  Specific  miRNAs  may  tend  to  cluster  together  suggesting  that  there  may  be 
patterns  of  miRNA  changes  that  need  to  be  explored  in  a  systems  approach; 

5)  Specific  alterations  in  oncogene  expression  or  loss  of  specific  suppressor  gene 
function  may  detennine  the  expression  pattern  of  miRNAs. 

Conclusions: 


The  results  from  this  study  provide  a  robust  dataset  of  genomic  information  at  three 
levels  -  miRNA,  mRNA,  and  genome  copy  number  -  from  mouse  models  of  mammary 
cancer  for  comparison  with  that  of  human  breast  cancer.  Analyses  of  the  data  indicate 
that  miRNA  patterns  are  quite  distinct  between  the  models  of  mammary  cancer 
dependent  upon  the  initiating  oncogenic  event.  This  may  in  part  be  due  to  the  cancer 
precursor  cell  that  is  being  targeted  by  the  manner  in  which  the  genetic  alteration  has 
been  engineered  in  the  mouse  germline,  or  it  may  be  due  in  part  to  the  particular  function 
of  the  oncogene  or  suppressor  gene  used  to  generate  the  tumors.  It  is  possible  that  part  of 
the  difference  in  miRNA  expression  between  the  models  is  therefore  detennined  by  the 
lineage  through  which  the  tumor  arises.  For  instance,  MMTV-her2/neu  and  MMTV-myc 
may  arise  through  the  targeting  of  tumorigenesis  to  precursor  cells  destined  to  become 
luminal  type  cells,  whereas  models  generated  through  the  disruption  of  p53  and/or 
BRCA1  function  may  target  a  more  primitive  precursor  cell.  Differences  in  cell  lineage 
could  also  account  for  changes  in  the  mRNA  expression  pattern. 

Ongoing  efforts  are  focused  on  further  identifying  the  major  similarities  in  miRNA 
expression  observed  in  the  mouse  mammary  tumors  with  that  reported  for  human  breast 
cancer.  miRNAs  meeting  this  criteria  will  be  further  studied  in  appropriate  in  vitro  and  in 
vivo  analyses. 


11 


REFERENCES 


1.  Iorio,  M.  V.,  Ferracin,  M.,  Liu,  C.  G.,  Veronese,  A.,  Spizzo,  R.,  Sabbioni,  S., 
Magri,  E.,  Pedriali,  M.,  Fabbri,  M.,  Campiglio,  M.,  Menard,  S.,  Palazzo,  J.  P.,  Rosenberg, 
A.,  Musiani,  P.,  Volinia,  S.,  Nenci,  I.,  Calin,  G.  A.,  Querzoli,  P.,  Negrini,  M.,  and  Croce, 
C.  M.  MicroRNA  gene  expression  deregulation  in  human  breast  cancer.  Cancer  Res,  65: 
7065-7070,  2005. 

2.  Alvarez-Garcia,  I.  and  Miska,  E.  A.  MicroRNA  functions  in  animal  development 
and  human  disease.  Development,  132:  4653-4662,  2005. 

3.  Lu,  J.,  Getz,  G.,  Miska,  E.  A.,  Alvarez-Saavedra,  E.,  Lamb,  J.,  Peck,  D.,  Sweet- 
Cordero,  A.,  Ebert,  B.  L.,  Mak,  R.  H.,  Ferrando,  A.  A.,  Downing,  J.  R.,  Jacks,  T., 

Horvitz,  H.  R.,  and  Golub,  T.  R.  MicroRNA  expression  profiles  classify  human  cancers. 
Nature,  435:  834-838,  2005. 

4.  He,  L.,  Thomson,  J.  M.,  Hemann,  M.  T.,  Hernando-Monge,  E.,  Mu,  D.,  Goodson, 
S.,  Powers,  S.,  Cordon-Cardo,  C.,  Lowe,  S.  W.,  Hannon,  G.  J.,  and  Hammond,  S.  M.  A 
microRNA  polycistron  as  a  potential  human  oncogene.  Nature,  435:  828-833,  2005. 

5.  Calin,  G.  A.,  Sevignani,  C.,  Dumitru,  C.  D.,  Hyslop,  T.,  Noch,  E.,  Yendamuri,  S., 
Shimizu,  M.,  Rattan,  S.,  Bullrich,  F.,  Negrini,  M.,  and  Croce,  C.  M.  Human  microRNA 
genes  are  frequently  located  at  fragile  sites  and  genomic  regions  involved  in  cancers. 

Proc  Natl  Acad  Sci  USA,  101:  2999-3004,  2004. 

6.  Deeb,  K.K.,  Michalowska,  A.M.,  Yoon,  C.,  Krummey,  S.M.,  Hoenerhoff,  M.J., 
Kavanaugh,  C.,  Li,  M-C.,  Demayo,  F.J.,  Linnoila,  I.,  Deng,  C-X.,  Lee,  E.Y-H.  P., 

Medina,  D.,  Shih,  J.H.,  and  Green,  J.E.  An  Integrated  Cancer  Genetic  Network  Predicts 
Aggressive  Human  Carcinomas  With  Poor  Prognosis,  Cancer  Research,  in  press. 

7.  Herschkowitz  JI,  Simin  K,  Weigman  VJ,  Mikaelian  I,  Usary  J,  Hu  Z,  Rasmussen 
KE,  Jones  LP,  Assefnia  S,  Chandrasekharan  S,  Backlund  MG,  Yin  Y,  Khramtsov  AI, 
Bastein  R,  Quackenbush  J,  Glazer  RI,  Brown  PH,  Green  JE,  Kopelovich  L,  Furth  PA, 
Palazzo  JP,  Olopade  OI,  Bernard  PS,  Churchill  GA,  Van  Dyke  T,  Perou  CM. 
Identification  of  conserved  gene  expression  features  between  murine  mammary 
carcinoma  models  and  human  breast  tumors. 

Genome  Biol.  2007  May  10;8(5):R76 


12 


Appendices 


13 


c-myc  Tu  04004022 _ whole  chr _ 

c-myc  Tu  04005648  3.4-82.1  whole  chr 


14 


Deletion  (Mb 


15 


16 


17 


chr  18  chr  X 

Amplification  (Mb)  I  Deletion  (Mb)  Amplification  (Mb)  I  Deletion  (Mb 


18 


